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In this paper, we show that two-dimensional billiards with point interactions inside exhibit a 
chaotic nature in the microscopic world, although their classical counterpart is non-chaotic. After 
deriving the transition matrix of the system by using the self-adjoint extension theory of functional 
analysis, we deduce the general condition for the appearance of chaos. The prediction is confirmed 
by numerically examining the statistical properties of energy spectrum of rectangular billiards with 
multiple point interactions inside. The dependence of the level statistics on the strength as well as 
the number of the scatterers is displayed. 

KEYWORDS: wave chaos, quantum mechanics, pseudointegrable billiard, point interaction, func- 
tional analysis 



I. INTRODUCTION 

The billiard is a suitable tool for examining generic fea- 
tures of dynamical systems. This is because, in spite of its 
simplicity, it covers the wide range of dynamical behav- 
iors, going from the most regular (integrable) to the most 
irregular (chaotic) depending on its shape. Pseudointe- 
grable |Q] billiards discussed in this paper is obtained 
by putting point scatterers into an integrable billiard. 
Their dynamical behavior is somewhat trivial in classical 
mechanics. A classical point particle never hits a point 
scatterer. More precisely, the set of trajectories of a point 
particle which hit point scatterers is of measure zero in 
the classical phase space. Thus we may neglect point 
obstacles in classical mechanics. However, the situation 
is drastically changed in quantum mechanics. Owing to 
the uncertain principle, the point particle gains a "size" 
and as a result, there might be a possibility of its hitting 
point scatterers. Our purpose is to show this is in fact 
the case and to clarify the condition for its occurrence. In 
order to emphasize that such phenomena originate from 
quantum effects, they might be called wave chaos ||. 

A particular case of a single scatterer located at the 
center of a rectangular billiard has been discussed in de- 
tails elsewhere ||- ||. We generalize the formulation 
and derive the transition matrix in case of multiple scat- 
terers. Based on this, a new numerical evidence for the 
existence of wave chaos is exhibited. We observe that the 
"degree" of chaos strongly depends on the number of the 
scatterers. 

The paper is organized as follows. In Sect. 2, we briefly 
mention the measure of chaos in quantum mechanics, 
which is described by statistics of energy levels. The 
quantum-mechanical treatment of point scatterers is dis- 
cussed in Sect. 3. In spite of seeming simplicity of the 
problem, a somewhat complicated argument based on 



self-adjoint extension is needed. In Sect. 4, we clarify the 
condition for the appearance of wave chaos and its va- 
lidity is confirmed by numerical experiments with a rect- 
angular billiard in Sect. 5. The paper is summarized in 
Sect.6. 



II. QUANTUM-MECHANICAL MEASURE OF 
CHAOS 

Although no mathematical proof exists, it is widely 
conjectured || that the statistical properties of quan- 
tum energy spectrum directly reflect the dynamical na- 
ture of the corresponding classical system. Generic in- 
tegrable systems obey Poisson statistics, while generic 
chaotic systems are described by the predictions of the 
Gaussian orthogonal ensembles (GOE) |Q. One of such 
statistical measures is the nearest-neighbor level spacing 
distribution P(S), which is defined such that P(S)dS 
gives the probability to find the spacing between any two 
neighboring energy eigenvalues in the interval (S, S+dS). 
Depending on the dynamical property of classical system, 
P(S) takes a form 



P(S) 



exp(-S), 



(integrable), 



ttS . irS 2 . . . . 
— exp( — ), (chaotic). 



(1) 



The degeneracy of eigenvalues frequently happens in inte- 
grable systems, while it hardly occurs in chaotic systems. 
Another convenient statistics is the spectral rigidity 



/ 1 f E+L/2 2 N 

A 3 (L) = ( min - / (N(E') ~ AE' - BY dE' 
\a,bLJ e _ l/2 ; 



(2) 



where N(E) is the staircase function which counts the 
number of eigenstates below energy E. The brackets (■) e 
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signify the average with respect to the energy in the en- 
ergy region under consideration. The As(L) statistics 
represents the average of the least square deviation of 
the staircase function N(E) from the best straight line 
fitting it in some intervals of length L. It takes a form 



A 3 (L) 



) 15' 



(integrable), 



r (lnL- 0.0687), (chaotic), 



(3) 



depending on the dynamical property. 



III. FORMULATION 

We first consider a quantum point particle of mass M 
moving freely in a two-dimensional bounded region (bil- 
liard) S. Let us denote the area of S by the same symbol. 
We impose the Dirichlct boundary condition so that wave 
functions vanish on the boundary of S. The eigenvalues 
and the corresponding normalized eigenfunctions are de- 
noted by E n and (p n (x) respectively; 



(4) 



The Hamiltonian Hq is the kinetic operator in L 2 (S) with 
domain D(H ) = H 2 (S) D H% (S) in terms of the Sobolev 
spaces. The Green's function of Ho is written as 



ip n (x)tp n {y) 



(•5) 



where to is the energy variable. The average level density 
of the system is given by 



MS 



(6) 



which is energy-independent. 

We now place a single point scatterer at x\ inside the 
billiard. The most naive manner for this purpose is to 
define the scatterer in terms of the Dirac's delta function; 



H = Ho + vS^ix-X!). 



(7) 



However, the Hamiltonian H is not mathematically 
sound. This can be seen from the eigenvalue equation 
of H, which is reduced to 



Vn{xi) 2 
UJ - E n 



(8) 



Since the level density is proportional to £?( d ~ 2 )/ 2 in spa- 
tial dimension d, the infinite series in Eq.(||) does not 
converge for d > 2. 



To handle the divergence, a scheme for renormalization 
is called for. The most general scheme is given by the self- 
adjoint extension theory of functional analysis ||]. We 
first consider in L 2 (S) the nonnegative operator 



H --*L 

Sl * 2M 



(9) 



with its closure H$ 1 in L 2 (S). Namely, we restrict D(H ) 
to the functions which vanish at the location of the point 
scatterer; 



D(HgJ = mx) G D(flo)|lK*i) = °}- 



(10) 



By using integration by parts, it is easy to prove that 
the operator Hg 1 is symmetric (Hermitian) . But it is not 
self-adjoint. Indeed, the equation for the adjoint of Hg 1 

H* Sl M*) = uil> u (2), ^ G D(H* Si ), (11) 

has a solution for Im ^ ||; 

^{x) = G (0 \x,x 1 -uu) 1 xeS-X!, (12) 

indicating 

D(H* Si ) = D(H Sl ) © N(H* Si -oj)® N(H* Si - Q) 

+ D(H Sl ), (13) 

where N(A) is the kernel of an operator A. Equation 
( |l3| ) means that H% 1 has the deficiency indices (1,1) and 
as a result, H$ 1 has one-parameter family of self-adjoint 
extensions He 1 ; 



D(H ei ) = {f\f = v + c(ip iA 



l ^-,a); 



^D(%),ceC,0<fl 1 <2i}, 



H 9l f = Hs t ip + iAc(ip 



iA 



(14) 



where A > is a scale mass. With the aid of Krein's 
formula ||, we can write down the Green's function for 
the Hamiltonian Hg 1 as 



G 6l (x,y;uj) = G^(x ) y^) 

+G<°> (x, X! ; lo)T 8i ( U )GM (X! , y; u) 



(15) 



In Eq. (jig) , the transition matrix Tg ± is calculated by 



where 



1 



(w - zA)c iA (w) - e ie i (w + iA)c_ lA (uj) ' 



c± iA (lu)= / G (0) (f,f 1 ;w)G (0) (f,f 1 ;±iA)dx. 
Js 

Using the resolvent equation, we have 
T Sl (a;) = («r 1 -G(a;))- 1 , 



(16) 



(17) 



(18) 
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where 



Wi -i = Acot^]T 



2 ^ El + A2 ' 

n— 1 n 



G{uj) = E Vn(xi) 2 ' 



E„ 



n=l 



y ui — E„ 



El + A 2 



(19) 
(20) 



The constant v\ is a coupling strength of the point scat- 
terer, the value of which ranges over the whole real num- 
ber as one varies < 6\ < 2n. It follows from Eq.(p"§|) 
that the eigenvalues of Hg ± are determined by 



G(w) 



(21) 



On any interval (E m , E m +\), the function G is monoton- 
ically decreasing, ranging over the whole real number. 
This means that the eigenvalue equation ( pl[ ) has a single 
solution cu m on each interval for any v\. The eigenfunc- 
tion of Hq 1 corresponding to an eigenvalue co m is given 
by 



ip m (x) <x G (0) (x, xi;uj ri 



oo , _ , 

<Pn{Xl) , . 

2^ 7. —fn( X )- 



n=l 



E n 



(22) 



We proceed to the case of multiple (N > 1) point scat- 
terers. Let us denote the position of the fc-th scatterer 
by x k . In this case, we should start with the operator 



H x = - 



2M 



(23) 



Cg°(S-X) 



and its closure Hx in L 2 (S). Here we set X = 
{x\, • ■ • , xn}. The operator Hx is symmetric (but not 
self-adjoint) and the equation 

H* x ^(x) = u^{x), i> u G D(H* X ), (24) 

has the N independent solutions for Im w^O; 

^\x) =G^{x,x k ;cu), xeS-X, (25) 

k = 1, • • • , N. This indicates that Hx has the deficiency 
indices (N, N) and as a result, it has 7V 2 -parameter fam- 
ily of self-adjoint extensions in general; 



iV 



N 



D(H UiX ) = {/|/ = <p + E c Ma + E CwCk)! 



fe=i 



<peD(H x ),c k eC}, 



N 



1=1 



N 



Hu,xf = H xV + i\ J2 Cfe ^iA " E U kiA), ( 26 ) 



k=l 



1 = 1 



where Uij denotes an arbitrary A^-dimcnsional unitary 
matrix. Using Krein's formula, we can relate the Green's 
function of H v x to G (0) ; 



Gu,x (x, y; u) = G (0) (x, y; w) 



N 



E G(° ) (f ) f^)r(.) M G<°)(x^», 



(27) 



k,l=l 



where the transition matrix T{u>) should satisfy 

- [Tfo/)]^ 
= G^(x k ,x l -Lu')~G {0 \x k ,x l -Lu). 



(28) 



Equation (£^) indicates that it is sufficient to define the 
transition matrix for some fixed u> since then T(u>) for any 
u> follows from Eq. (p8|) . Also note that Eq.(pTj) implies 



T( W )t=T(o). 

The general form of T{w) which satisfies Eqs. 
(p9|) is given by 



Tiujy 1 = A-G(u), 



(29) 
and 

(30) 



where 



oo _^ 

E fn(x k ) 2 ( 

G(Uj) k l = < «=! 



E n 



LO - E n 



El+A 2 ' 
(k = l), 



(31) 



and A is any A^-dimensional Hermitian matrix. If A is 
not diagonal, different scatterers x k are connected by the 
boundary conditions. On the other hand, if A is diag- 
onal the point scatterers are independent, namely, one 
has separated boundary conditions at each x k . We here 
restrict ourselves to the latter which is important from a 
practical point of view. In this case, A can be parame- 
terized as 



A kl = 5, 



kiv k 



v k G R. 



(32) 



One can regard v k as the strength of the fc-th scatterer. 
The eigenvalue of Hu x is determined by the poles of 

det(T(u)- x ) = 0, (33) 

which is reduced to Eq.(pT|) for N = 1. 

A final remark is that one can show the relation 

[/ = - i [T(-iA)T(iA)" 1 ]. (34) 

Equation (Q) indicates that the unitarity of U is equiv- 
alent to the fact that T{ih) is a normal matrix. 



IV. CONDITION FOR STRONG COUPLING 

In this section, we discuss the condition for strong cou- 
pling under which point scatterers have a substantial ef- 
fect on energy spectrum of the empty billiard. For this 
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purpose, we examine the energy-dependence of the inflec- 
tion points of G in Eq. ( |2~C| ) , based on a semi-quantitative 
argument suitable for examining the statistical properties 
of spectrum. 

We first consider the case of a single scatterer. The 
first notice is that the average value of ip n (xi) 2 among 
many n is independent of the energy; 



(<Pr 



1/5. 



(35) 



2|) that if Lj m ~ E m (resp. 
ip m (resp. tp 



We thus recognize from Eq 
E m+ i) for some m, then 
This implies that a point scatterer distorts the wave func- 
tion if the eigenvalue to m is located around the midpoint 
of the interval (E mi E m+ i). For such u) m , the value of 
G can be estimated by the principal integral as follows, 
since the contributions on the summation of G from the 
terms with n ~ m cancel each other; 



G(uj m ) ~ g{u m ), uj„ 



g(uj) = (ip n (x 1 ) 2 ) n p av P 

__ M u) 
_ 27r n A' 



E m + E, 



m+l 



(36) 



1 



- E E 2 + A 2 



dE 



(37) 



Here we have defined a continuous function g(u>) which 
behaves like an interpolation of the inflection points of G. 
Equation (|36| ) indicates that the wave function mixing 
mainly occurs in the eigenstate with an eigenvalue such 
that 



/(w m ) 



(38) 



The "width" of the strong coupling region (allowable 
error of v^ 1 in Eq.(|38|)) is estimated by considering a 
linearized eigenvalue equation. Expanding G at oj m ~ 
(E m + E rn+ i)/2, we can rewrite the eigenvalue equation 
as 



{Vn{xi) 2 ) n 



1 



-8(^„(f 1 ) 2 ) n pLE ( 2n -l) 2 

n=l y ' 

= w 2 (<p n (Zi) 2 ) n plv (42) 

The second equality follows from the approximation that 
the unperturbed eigenvalues are distributed with a mean 
interval in the whole energy region. Noticing Eqs.(|^) 
and (|35|), we obtain 



7fM 

~2~' 



(43) 



which is independent of the energy lo. We can summarize 
the findings as follows; The effect of a point scatterer of 
coupling strength v\ is substantial mainly in the eigen- 
states with eigenvalue to such that 



M , LU 
— In — 

2tt A 



7rM 



(44) 



The condition ( J44| ) is generalized to the case of multiple 
scatterers. Equation ([55]), which is the main assumption 
in the above argument, is valid even for the perturbed 
eigenfunction (p2f). This allows us to repeat the above 
when the number of the point scatterers increases. We 
conclude that the fc-th point scatterer with strength 
affects the energy spectrum in the energy region such 
that 



M , u 
— In — 
2tt A 



, ttM 

$ , k = l,---,N. 

~ 4 ' ' ' 



(45) 



It is worthy to emphasize that the strong coupling re- 
gion shifts with a logarithmic dependence of energy. This 
indicates that the system recovers the integrability in the 
semiclassical (high-energy) limit for any strength Vk- 



G'(u) m )(uj - u m ) ~ v 1 -G(u m ). 



(39) 



V. NUMERICAL EXAMPLE 



In order that Eq. (|39|) has a solution u> ~ w m , the range 
of RHS has to be restricted to 



G(ui ri 



< * 



A=\G'(uj m )\p : 



-1 

av i 



(40) 



(41) 



where we have defined the width A which is nothing 
but the variance of the linearized G on the interval 
(E m , E m+ i). The magnitude of \G'(uo m )\ can be esti- 
mated as follows; 



TABLE I. Location of the point scatterers 



k 




Vk 


1 


0.6224826 


0.2758356 


2 


0.8202505 


0.4561782 


3 


0.1802603 


0.6365209 


4 


0.3780281 


0.8168635 


5 


0.5757960 


0.2332624 


6 


0.7735638 


0.4136051 


7 


0.1335736 


0.5939477 


8 


0.3313415 


0.7742904 


9 


0.5291093 


0.1906893 


10 


0.7268772 


0.3710320 
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In this section, we set the scale mass A = 1 with- 
out losing generality. This makes all physical quanti- 
ties dimensionless. We consider a quantum particle with 
mass M = 2~k moving in a rectangular billiard with side- 
lengths (l x ,l y ) — (7r/3, 3/tt), and hence 5 = 1. The av- 
erage level density is p av = 1, according to Eq.(^). The 
eigenvalue E nsc ^ ly and the corresponding eigcnfunction 
^n x ,n„ are given by 




P(S) 0.6 - 



P(S) 0.6 j 



0.5 1 1.5 2 2.5 3 

s 

FIG. 1. The dependence of the nearest-neighbor level spac- 
ing distribution on the strength of the point scatterers. The 
number of the scatterers is N = 5 in all cases. The strength 
of all the scatterers is common. The statistics is taken among 
the eigenstates between uioo and Wiioo- The solid (broken) 
line is the prediction of GOE (Poisson statistics). 



E = 




+ 



4 . n x irx . n v iry 
sm— — sin^ — , 



(46) 
(47) 



respectively, (n x ,n y — 1,2, •••). We obtain E n and <p n 
by arranging the double-indexed eigenvalues and eigen- 
functions in ascending order of magnitude of E nx ^ n . 



A 3 (L) 




2 4 6 8 10 12 14 16 18 20 

L 

FIG. 2. The dependence of the spectral rigidity on the 
strength of the point scatterers. The other indications are 
the same as in Fig.l. 



A 3 (L) 




2 4 6 8 10 12 14 16 18 20 

L 

FIG. 3. The dependence of the spectral rigidity on the 
number of the point scatterers. The strength of all the scat- 



terers is v. 



7.5. 



We put several point scatterers inside the rectangle. 
The location of point scatterers, Xk — {%k,yk), is shown 
in Table The first N coordinates are used in case of 
the number of scatterers being N. It is expected from 
Eq.([45"|) that the fc-th scatterer has a substantial effect 
on energy spectrum if its strength satisfies 



In a; < 



7T" 

T 



4.9348. 



(48) 



This is in fact the case, as seen from Fig.|] where the 
nearest- neighbor level spacing distribution P(S) for three 
values of the strength of the scatterers is shown. The 
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number of point scatterers is N = 5 in all cases. The 
strength of all the scatterers is common. The statis- 
tics is taken within one thousand eigenvalues between 
ujioo = 110.3579 and Wxioo = 1138.9682. In this energy 
region, the chaotic spectrum is expected to appear for 
v^ 1 ~ 5 ~ 7. The coincidence with the GOE prediction 
(solid line) is satisfactory for v^ 1 — 5. As v^ 1 increases, 
the level repulsion becomes weak and P(S) shows inter- 
mediate between Poisson and GOE for v^ 1 = 12.5, as 
expected from Eq.(^8|). A further increase of v^ 1 makes 
the distribution approach the Poisson distribution (bro- 
ken line), though there still remains the level repulsion 
even for v^ 1 = 20. A similar tendency can be observed in 
the spectral rigidity A3 (£) , which FigJ| shows for several 
values of v^ 1 under the same condition as in Fig.0. The 
A^(L) statistics is close to the GOE prediction (solid 
line) in case of v^ 1 — 5 or 7.5. As v^ 1 increases, the 
gradual approach to the Poisson statistics (broken line) 
is observed. There still remains a considerable difference 
from the GOE prediction even in case of the maximal 
coupling. However, it tends to disappear as the number 
of scatterers increases. Figure |^ shows the dependence 
of the spectral rigidity on the number of the scatterers. 
The strength of all the scatterers is v^ 1 = 7.5, which is 
close to the maximal coupling strength. One can observe 
the gradual shift to the GOE prediction as the number 
of the scatterers increases. 

VI. CONCLUSION 

We have discussed the condition for the appearance of 
wave chaos in the two-dimensional quantum pseudoin- 
tegrable billiards with multiple point scatterers inside. 
Chaotic spectrum appears if the condition of Eq. ( Ma) is 
satisfied. It is described by a logarithmically energy- 
dependent strip with a constant width. The validity 
of our prediction is confirmed by numerical experiments 
with a rectangular billiard which contains multiple point 
obstacles inside. The degree of chaos depends on the 
number of the scatterers. We observe the GOE-like spec- 
trum for N = 10 scatterers with the maximal coupling 
strength. 

It should be stressed that, besides a fundamental as- 
pect as a dynamical system, the quantum billiard with 
point interactions has a close relation to real systems. 
The quantum billiard is a natural starting point for ex- 
amining the particle motion in microscopic bounded re- 
gions. The rapid progress in the microscopic or meso- 
scopic technology makes it possible to realize such set- 
tings. Real systems are, however, not free from small 
impurities which affect the particle motion inside. In 
the presence of a small amount of contamination, even a 
single-electron problem becomes unmanageable in an an- 
alytic manner. The modeling of the impurities with point 
interactions makes the problem easy to handle without 



changing essential dynamics. We hope that the current 
work serves as a guideline of the future research on the 
effect of impurities in mesoscopic devices. 



[1] P.J. Richens and M.V. Berry, "Pseudointegrable systems 
in classical and quantum mechanics," Physica D, vol.2, 
pp.495-512, 1981. 

[2] S. Albeverio and P. Seba, "Wave chaos in quantum 
systems with point interaction," J. Stat. Phys., vol.64, 
nos.1/2, pp. 369-383, 1991. 

[3] T. Shigehara, "Conditions for the appearance of wave 
chaos in quantum singular systems with a pointlike scat- 
terer," Phys. Rev. E, vol.50, no.6, pp.4357-4370, Decem- 
ber, 1994. 

[4] T. Shigehara and T. Cheon, "Spectral properties of three- 
dimensional quantum billiards with a pointlike scatterer," 
Phys. Rev. E, vol.55, no.6, pp.6832-6844, June, 1997. 

[5] T. Shigehara, H. Mizoguchi, T. Mishima, and T. Cheon, 
"Wave chaos in quantum pseudointegrable billiards," 
Proc. International Symposium on Nonlinear Theory and 
its Applications, NOLTA'97, pp.441-444, Honolulu, Nov., 
1997. 

[6] O. Bohigas, "Random matrix theories and chaotic dy- 
namics," in Les Houches Summer School Proceedings Ses- 
sion 52, ed. M.-J. Giannoni, A. Voros, and J. Zinn- Justin, 
pp.87-200, North-Holland, Amsterdam, 1989. 

[7] M.L. Mehta, "Random matrices and the statistical theory 
of energy levels," Academic Press, New York, 1967, revised 
and enlarged edition, 1990. 

[8] S. Albeverio, F. Gesztesy, R. H0egh-Krohn, and H. 
Holden, "Solvable models in quantum mechanics," 
Springer- Verlag, New York, 1988. 

[9] J. Zorbas, "Perturbation of self-adjoint operators by Dirac 
distributions," J. Math. Phys., vol.21, no.4, pp.840-847, 
April, 1980. 



G 



Profile of the Authors 



Takaomi Shigehara 

email : sigehara@ics . sait ama-u . ac . jp 

http: / /www. me. ics.saitama-u.ac.jp/~sigehara/ 
received the B.S, M.S and Ph.D. degrees in Physics from 
the University of Tokyo in 1983, 1985 and 1988, respec- 
tively. He is currently an Assistant Professor in the 
Department of Information and Computer Sciences at 
Saitama University. His research interests are quantum 
chaos, high-performance computing, and numerical anal- 
ysis. 



Hiroshi Mizoguchi 

email: hm@ics.saitama-u.ac.jp 

http: / /www. me. ics.saitama-u.ac.jp/~hm/ 

received the B.E. degree in Mathematical Engineering in 
1980, and the M.E. and Ph.D. degrees in Information 
Engineering in 1982 and 1985, respectively, from the Uni- 
versity of Tokyo. He is currently an Associate Professor 
in the Department of Information and Computer Sci- 
ences at Saitama University. His research interests are 
vision processor system, robotics, and quantum chaos. 



Takctoshi Mishima 
email: misima@ics.saitama-u.ac.jp 
http: / /www. me. ics.saitama-u.ac.jp/~mishima/ 
received the B.E., M.E. and Ph.D. degrees in Electrical 
Engineering from Meiji University, Japan, in 1968, 1970 
and 1973, respectively. He is currently a Professor in 
the Department of Information and Computer Sciences 
at Saitama University. His research interests are foun- 
dation of symbolic and algebraic computation, axiomatic 
logic system, mathematical pattern recognition, and quan- 
tum chaos. 



Taksu Cheon 

email: cheonQ mech . ko chi- 1 ech . ac . j p 

http: / /www. kochi-tech.ac.jp/~cheon/ 
received the B.S., M.S. and Ph.D. degrees in Physics 
from the University of Tokyo in 1980, 1982 and 1985, 
respectively. He is currently an Associate Professor in 
the Laboratory of Physics at Kochi University of Tech- 
nology. His research interests are quantum mechanics, 
chaos, and quantum chaos. 



7 



